1
Dasar Teori dan Kelayakan Masalah
MATH007Lesson 5
00:00
Sebelum kita menggunakan kekuatan solver numerik seperti metode Runge-Kutta orde keempat atau formula Adams-Moulton, kita harus menanyakan pertanyaan mendasar: Apakah solusi benar-benar ada, dan apakah stabil? Dasar teori masalah nilai awal (IVP) memberikan lampu hijau matematis, memastikan bahwa diskretisasi kita konvergen menuju realitas fisik yang bermakna daripada kebisingan numerik.

Dasar Utama: Kontinuitas Lipschitz

Untuk mengendalikan bagaimana kesalahan menyebar, kita membutuhkan fungsi $f(t, y)$ yang tidak "melompat" terlalu liar. Hal ini diformalisasikan oleh Kondisi Lipschitz.

Definisi 5.1: Kondisi Lipschitz

Suatu fungsi $f(t, y)$ memenuhi kondisi Lipschitz dalam variabel $y$ pada himpunan $D \subset \mathbb{R}^2$ jika suatu konstanta $L > 0$ ada dengan:

$$|f(t, y_1) - f(t, y_2)| \le L|y_1 - y_2|$$

untuk semua $(t, y_1), (t, y_2) \in D$. Konstanta ini $L$ adalah "batas kecepatan" untuk perubahan vertikal fungsi tersebut.

Contoh 1: Analisis Konstanta Lipschitz

Pertimbangkan $f(t, y) = t|y|$ pada $D = \{(t, y) \mid 1 \le t \le 2, -3 \le y \le 4\}$. Berdasarkan Teorema Nilai Rata-Rata (atau sifat nilai absolut):

$|f(t, y_1) - f(t, y_2)| = |t|y_1| - t|y_2|| = |t| \cdot ||y_1| - |y_2|| \le |t| \cdot |y_1 - y_2|$.

Karena nilai maksimum $t$ di domain kita adalah 2, maka konstanta Lipschitznya adalah $L=2$.

Integritas Geometris Domain

Kita tidak dapat menyelesaikan IVP di domain yang penuh lubang. Kita membutuhkan Kekonveksan.

Definisi 5.2: Himpunan Konveks

Suatu himpunan $D$ dikatakan konveks jika untuk dua titik sembarang $(t_1, y_1)$ dan $(t_2, y_2)$, segmen yang didefinisikan oleh:

$$((1 - \lambda)t_1 + \lambda t_2, (1 - \lambda)y_1 + \lambda y_2)$$

untuk $\lambda \in [0, 1]$ juga terkandung dalam $D$. Ini memastikan tidak ada bagian dari lintasan solusi "keluar" dari zona komputasi yang valid.

Teorema Kehadiran dan Kepemilikan Unik

Ketika kondisi-kondisi ini sesuai, kita menggali Teorema 5.4: Jika $f$ kontinu dan memenuhi kondisi Lipschitz pada himpunan konveks $D$, maka IVP $y' = f(t, y), y(a) = \alpha$ memiliki unik solusi $y(t)$. Ini menjelaskan metode se-sederhana Euler ($w_{i+1} = w_i + h f(t_i, w_i)$) atau se-sekompleks logika prediktor-korektor:

$WC = w_{i-1} + \frac{h}{24}[9f(t_i, WP) + 19f(t_{i-1}, w_{i-1}) - 5f(t_{i-2}, w_{i-2}) + f(t_{i-3}, w_{i-3})]$.

🎯 Prinsip Utama: Kelayakan Masalah
Suatu masalah dikatakan layak jika solusi unik ada dan bergantung secara kontinu pada data awal. Jika konstanta Lipschitz $L$ sangat besar, masalah menjadi "keras". Pada persamaan keras, bagian transien cepat meluruh, tetapi turunannya (magnitudo $c^n e^{-ct}$) tidak, sehingga diperlukan Algoritma 5.8: Trapesium dengan Iterasi Newton untuk menjaga stabilitas.